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We envision that quantum spin Hall effect should be observed in (lll)-oriented thin 
films of SnSe and SnTe topological crystalline insulators. Using a tight-binding ap¬ 
proach supported by first-principles calculations of the band structures we demonstrate 
that in these films the energy gaps in the two-dimensional band spectrum depend in 
an oscillatory fashion on the layer thickness. These results as well as the calculated 
topological invariant indexes and edge state spin polarizations show that for films ~ 
20—40 monolayers thick a two-dimensional topological insulator phase appears. In this 
range of thicknesses in both, SnSe and SnTe, (lll)-oriented films edge states with 
Dirac cones with opposite spin polarization in their two branches are obtained. While 
in the SnTe layers a single Dirac cone appears at the projection of the P point of the 
two-dimensional Brillouin zone, in the SnSe (lll)-oriented layers three Dirac cones at 
M points projections are predicted. 

PACS numbers: 71.20.-b, 71.70.Ej, 73.20.At, 79.60.-i 


Topological insulators (TIs) and the quantum spin Hall 
(QSH) effect attract significant interest for both funda¬ 
mental and practical reasons. In the three-dimensional 
(3D) and two-dimensional (2D) TIs the bulk insulating 
states are accompanied by metallic helical Dirac-like elec¬ 
tronic states on the surface (edges) of the crystal. Due 
to the time-reversal symmetry, these surface states are 
topologically protected against scattering at 180°. In 
2D TIs this means that the metallic edge states pro¬ 
vide dissipationless, spin-polarized conduction channels. 
Such property makes these structures extremely inter¬ 
esting for low-power-consumption electronics and spin- 
tronics. Thus, the quest for systems with topologically 
non-trivial edge states, which give rise to QSH effect, has 
become recently one of the most important topics in con¬ 
densed matter physics. 

The QSH effect is the net result of two opposite polar¬ 
ized spin currents traveling in opposite directions along 
the edges of a 2D TI. It has been predicted by Kane and 
Mele in graphene. 1 In this material, however, it could 
occur only at unrealistically low temperatures, since the 
intrinsic spin-orbit coupling, which should open a band 
gap at the Dirac points, is very weak. Subsequently, 
it has been proposed and confirmed experimentally that 
QSH state might arise in HgTe/HgCdTe quantum wells 2 
and also in InAs/GaSb heterostructures. 3 In these ex¬ 
periments the transition from insulating to conducting 
behavior has been observed. In both cases the conduc¬ 
tivity was close to the 2 e 2 /h value expected for the two 
parallel quantum Hall channels. The spin polarized na¬ 
ture of the edge states has been shown much later and 
only in the HgTe-based structure. 4 All these observations 
were made at very low temperatures, below 10 K, due to 
the small energy gap in these systems. 

In the search for new QSH structures many theoretical 
predictions have been made, in which various, sometimes 


exotic, chemical classes of 2D materials were considered. 
These include slightly buckled honeycomb lattice of Si 
atoms (silicene), 5 Bi bilayers on different substrates, 6 
Ge(BL c Sbi_ a; ) 2 Te 4 with various Bi concentrations 7 and 
functionalized ultrathin tin films (stanene). 8 Recently, it 
has been also shown that 2D transition metal dichalco- 
genides (Mo or W) should form a new class of large-gap 
QSH insulators. 9 Finally, a 2D TI with energy gap as 
large as 0.8 eV has been predicted for an overlayer of Bi 
grown on semiconductor Si(lll) surface functionalized 
with one-third monolayer of halogen atoms. 10 Despite all 
these theoretical predictions, no stand-alone thin film or 
a thin film supported on a suitable substrate have been 
up to now experimentally demonstrated to harbor a QSH 
state. 

In this article we describe a study of (lll)-oriented 
thin films of well known IV-VI semiconductors SnSe and 
SnTe, which can be quite easily grown on a BaF 2 sub¬ 
strate by e.g. molecular beam epitaxy. Such layers 
have rock-salt crystal structure, despite that bulk SnSe 
has an orthorhombic structure. 11 Recently, it has been 
shown that these compounds are so-called topological 
crystalline insulators (TCIs). TCIs are nontrivial insu¬ 
lators supporting surface Dirac fermions protected not 
by time-reversal but by crystal symmetry. 12 The IV-VI 
semiconductors, i.e., SnSe and SnTe, as well as PbTe and 
PbSe and the substitutional solid solutions of Pb- and Sn- 
based chalcogenides were identified earlier as trivial insu¬ 
lators. This is because in these compounds the band in¬ 
version happens simultaneously at an even number (four) 
of L points in the Brillouin zone (BZ). The angle-resolved 
photoemission spectroscopy (ARPES) confirmed, how¬ 
ever, that these materials belong to the TCI class. It 
has been shown that metallic surface states exist on the 
(001) surfaces of SnTe 13 as well as Pbi-xSn^Se 14 and 
Pbi-^Sn^Te. 15 These gapless (001) surface states are 
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supported by mirror symmetry. It should be empha¬ 
sized that due to the spin-orbit interactions the metal¬ 
lic surface states observed in real SnTe- or SnSe-based 
compounds have almost linear, Dirac-like dispersions. It 
has been shown lately that in (OOl)-oriented thin films of 
TCIs hybridization leads to a new 2D TCI phase, which 
supports two pairs of spin-filtered edge states. 16,17 These 
edge states are protected solely by mirror symmetry - the 
topological phase is indexed by the mirror Chern number 
\N m \ =2. In contrast, our calculations show that in ul- 
trathin TCI films (SnTe and SnSe) grown not along (001) 
but along (111) crystallographic axis, the QSH state can 
be obtained. This is because the four L point projec¬ 
tions onto the 2D BZ of a (lll)-oriented thin film are not 
equivalent. While the energy structures at the three M 
projections are the same by symmetry, the energy struc¬ 
ture at the fourth projection at the r point is different. 
This results in a possibility that the band inversion takes 
place in an odd (either one or three) number of points. 
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FIG. 1: Energy gaps in (lll)-oriented films with 
an odd number of monolayers. The dependence of 
the 2D states band gaps on the thickness of the (a) 
SnTe and (b) SnSe cation-terminated film. The green 
shadowed areas denote thicknesses, for which the value 
v = 1 has been obtained. 



A IV-VI crystal in rock-salt structure has {111} planes 
alternately composed of either cations or anions. Thus, 


(lll)-oriented slabs consisting of an even number of 
monolayers have one surface cation- and the other anion- 
terminated. In contrast, in slabs with an odd number of 
layers both surfaces are the same and the inversion sym¬ 
metry is preserved. These two cases have to be distin¬ 
guished. 


I. ODD NUMBER OF MONOLAYERS 


For a (lll)-oriented SnTe or SnSe thick film four sin¬ 
gle, topologically protected Dirac-cones in the four pro¬ 
jections of the L-points onto the 2D BZ (one in V point 
and three in the M points) are obtained in the calcula¬ 
tions. This is similar to what we have obtained before for 
a (lll)-oriented bulk PbSnTe crystal. 18 For the anion- 
terminated surfaces, the bands are brought to contact 
forming anion Dirac cones, while in the other case the 
bands meet to form cation Dirac cones. As shown already 
in Ref. 18 and 19, all L points belong to the three {110} 
mirror planes of the (lll)-surface and all Dirac points 
should be thus topologically protected. This predictions 
have been also confirmed experimentally for (Pb,Sn)Se 20 
as well as for SnTe 21 case. While at the r point the 
Dirac-cone is isotropic, the band structure around M is 
strongly anisotropic, i.e., depends differently on the k val¬ 
ues along the M - T and M - K directions. This is due 
to different orientations of the constant energy ellipsoids 
around different L points: the L i ellipsoid is projected on 
r along its long axis, whereas the ellipsoids for the other 
three L points are inclined away from the projection di¬ 
rection. When the thickness of such TCI film decreases 
the gaps in the Dirac cones open, since the wave functions 
of the top and bottom surface states start to hybridize 
with each other. Examples of band structures around 
r and around M points in the 2D BZ of (lll)-oriented 
SnTe(SnSe) thin films are presented in Fig. SI (Fig. S2) 
in the Supplement. In Supplementary Section I the cal¬ 
culation method together with the obtained tight-binding 
parameters is also presented. Here in Fig. 1 we present 
the 2D state energy gaps in ultrathin cation-terminated 
(lll)-oriented SnTe (a) and SnSe (b) slabs as a function 
of the number (odd) of monolayers in the film. 

As shown in Fig. 1, in SnTe films the energy gap at 
the M point decreases monotonically with the increase 
of the film thickness. At the r point, however, damped 
oscillations of the energy gap with the number of mono- 
layers can be observed. According to these results the 
first and biggest band gap inversion can be expected for 
SnTe films ~ 3 - 6 nm thick, i.e., consisting of 17-33 
monolayers. In contrast, in SnSe thin films energy gap 
at r decreases monotonically, while the one at the M 
point oscillates with the thickness of the film - in this 
case the first band gap inversion in M should occur for 
films with 19-39 monolayers, what corresponds to sim¬ 
ilar thicknesses of 3.3 - 6.7 nm. Analogous results are 
obtained for the anion-terminated layers. In this case, 
however, due to the smaller range of anion wave func- 
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tions the inverted band gap starts to occur for slightly 
thinner layers (2.4 - 5.3 nm for SnTe and 2.6 - 5.7 nm for 
SnSe), as shown in the Supplementary Section II. The 
damped oscillations of the band gap result from the k 2 
corrections to the Dirac-like Hamiltonian of massive elec¬ 
tron and have been already predicted for thin films of 
3D TIs, like Bi 2 Se 3 or Bi 2 Te 3. 22,23 In SnTe and SnSe k 2 
terms are highly anisotropic. The differences in the re¬ 
sults obtained for films of these materials stem from the 
fact that in SnTe the k 2 term moves the energy gap from 
L towards the r point, whereas in SnSe towards the W 
point (see Fig. S3 in the Supplement). 

In analogy to results obtained in Ref. 17 for (001)- 
oriented Pbi-^Sn^Te layers, the described above changes 
of the energy gaps are associated with a topological phase 
transition. Our calculations show, however, that in the 
case of (lll)-oriented SnTe and SnSe layers the topo¬ 
logical phase changes from trivial insulator not to TCI, 
but to the 2D TI. The effect is similar to the transition 
to 2D TI predicted for Bi 2 Se 3 or Bi 2 Te 3 thin films . 23 
This result is obtained by calculations of the invari¬ 
ant v = 0,1, which in 2D distinguishes between the trivial 
and QSH phase. For systems with inversion symmetry, 
i.e., for films consisting of an odd number of monolayers, 
the parity of the occupied bands eigenstates at the four 
time-reversal momenta points in the BZ (one F and three 
M) determines the Zi invariants . 24 The green shadowed 
areas in Fig. 1 denote the SnTe and SnSe layer thick¬ 
nesses, for which the value v = 1 has been obtained by 
the prescription given in Ref. 24. 

In Fig. 2a and 2b we present the calculated [110] edge 
spectral functions of the (lll)-SnTe cation-terminated 
film consisting of 19 monolayers (i.e., with inverted band 
gap - compare Fig. la). The results for a slab consisting 
of 13 monolayers, where a trivial insulator phase is ex¬ 
pected, is presented in Fig. 2c. The 2D BZ zone of the 
(lll)-oriented film and ID BZ of the [110] edge are pre¬ 
sented in Fig. 2d. As shown in the latter, the F and Mi 
points project onto the edge at k x = 0, while M 2 and 
M 3 at k x = iy^Tr/ao- In agreement with the results 
of topological invariant calculations, shown above, Dirac 
crossing of the edge states appears in the band gap at the 
r projection for the 19-layers, but not for the 13-layers 
thick slab (compare Fig. 2a and 2c). At the projection 
of M 2 and M 3 (for k x = y/Onr/ao) the Dirac cones are 
not obtained for any SnTe layer thickness. In contrast, in 
SnSe, as shown in Fig. 3a and 3b for the anion-terminated 
21-monolayer thick slab, two Dirac crossings of the edge 
states appear inside the band gap at the k x = \Z2i:/ao 
point. As the third Mi point projects at k x = 0, another 
edge state with a Dirac node appears at k x = 0. Due 
to the strong band overlapping, the latter is dispersed 
within the valence band. 

Finally, we have calculated the G-j- and Gj, contribu¬ 
tions of the spin up and spin down (lll)-projections in 
the edge spectral functions. The sign of the G-f - Gj. 
difference corresponds to spin polarization of the edge 
states. As expected for the 2D TI in the QSH state, 



FIG. 2: Edge states in cation-terminated 
(lll)-SnTe film. The calculated spectral functions for 
a slab of 19 monolayers at the projection of (a) r and 
(b) M 2 and M 3 points onto the [110] edge, (c) The 
spectral function at the projection of r for 13 
monolayers thick film, (d) The 2D BZ of (lll)-oriented 
film (green hexagon) and ID BZ of its [110] edge (red 
line). Local extrema of the 2D bands projected to the 
edge are here and in the following figures denoted by 
black dotted lines. 


the obtained edge states are spin polarized and the spin 
polarizations in the two branches of the Dirac cones are 
always opposite. These spin polarizations in all studied 
cases are about 80%, i.e., the contribution of the domi¬ 
nant spin configuration surpasses ca 8-10 times the op¬ 
posite spin orientation. 


II. EVEN NUMBER OF MONOLAYERS 

Unlike the case of an odd number of layers, for (111)- 
oriented slabs consisting of an even number of monolay¬ 
ers there is no inversion symmetry in the system and the 
energy gaps are moved from the time-reversal invariant 
points (see Figs Sl(b) and S2(b) in the Supplementary 
Section I). The gap inversion occurs now in two points 
related by the time reversal symmetry on two sides of 
each time reversal invariant k point. It has been shown 
for inversion-asymmetric systems in general 25 that the 
gap closing can be accompanied by transition between 
the trivial insulator and the QSH phase when the gap 
inversion occurs in an odd number of such pairs. In 
SnSe films this takes place in three pairs in the vicin¬ 
ity of the three nonequivalent M points while in SnTe 
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FIG. 3: Edge states in 21-monolayers thick 
anion-terminated (lll)-SnSe film. The edge 
spectral functions in the vicinity of (a) k x = 0 (a) and 
(b) k x = y/2ir/ao points of the ID BZ. The 
corresponding spin densities are shown on the right. 
Red and green lines represent the spin-down and 
spin-up polarization, respectively. 


in one pair around the F point. Our calculations show 
that indeed for a given range of thicknesses the topolog¬ 
ical phase transition to the TI phase can be expected in 
these (lll)-oriented films. 

If there is no inversion symmetry in the system, the 
Z 2 topological invariants can not be determined by the 
method utilizing the parity of the eigenfunctions, which 
we used for films with an odd number of monolayers. To 
calculate Z 2 for such layers we have adopted the method 
proposed by Fukai and Hatsugai. 26 In this method the 
improvements to computing Chern numbers in a lattice 
BZ enable the tight-binding calculations of Z 2 based on 
the formula derived by Fu and Kane in Ref. 27. The de¬ 
pendence of the band gap of SnTe (lll)-oriented films 
consisting of an even number of monolayers on the layer 
thickness is shown in Fig. 4a. The thicknesses for which 
the value v = 1 has been obtained are marked by the 
green shadow. As one can see in Fig. 4a the QSH phase 
can be expected in SnTe films consisting of 14-28 mono- 
layers (~ 2.5 - 5.1 nm thick). For SnSe the QSH phase 
is obtained for 18-38 monolayers, i.e., for film thicknesses 
in the range 3.1 - 6.6 nm. In the case of SnTe films with 
even number of layers the whole inverted energy gap for 
states in the vicinity of F is situated inside the large gap 
between conduction and valence bands close to the M 
point, as one can clearly see in Fig. 4b. The largest value 



FIG. 4: Energy gaps in SnTe film with an even 
number of monolayers, (a) The dependence of the 
2D state band gaps on the slab thickness. For 
thicknesses within the green shadowed area v = 1 was 
obtained, (b) The valence and conduction bands 
extrema in the vicinity of F (black) and M (red) vs the 
number of monolayers. 


of the inverted gap, about 75 meV, is obtained for the 
18-20 monolayers thick SnTe slab. 

The calculated spectral functions along the [110] edge 
for the (lll)-oriented, 18-monolayers thick SnTe film 
show a clear Dirac node in the center of the band gap 
at the k x = 0 point. This is shown in Fig. 5a. The ob¬ 
tained results are in full agreement with the formalism 
developed in Ref. 25. We note that a similar result has 
been obtained for graphene 1 - in this case one observes 
a single edge state with a Dirac node at F (time-reversal 
symmetry point), despite the energy gaps appear in two 
points K (in which the time-reversal is not preserved). 

Fig. 6 presents the edge spectral functions at the pro¬ 
jections of M points at (a) k x = 0 and (b) k x = \/2Tr/ao 
point, calculated for the 20-monolayers thick SnSe slab. 
As shown in the Supplementary Material, in this case the 
band gap at the k x = 0 is very small. Still, a pair of edge 
states can be seen near the top of the valence band. The 
spin polarization of this edge state is shown in Fig. 6c. At 
the k x = \/2ir/a$ point, however, the edge states which 
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FIG. 5: Edge states in 18 monolayers thick SnTe 
film, (a) The edge spectral functions at the vicinity of 
k x = 0 point of the ID BZ. (b) The corresponding spin 
density. The spin-down polarization is denoted by red, 
the spin-up by green color. 




FIG. 6: Edge states in 20-monolayer thick SnSe 

film. The edge spectral functions at k x = 0 are shown 
in (a), at the k x = \/2ir/a$ point in (b). The spin 
polarizations of the edge states presented in (a) and (b) 
are shown in (c) and (d), respectively. Again, green 
color represents spin-up and red spin-down contribution 
to the edge spectral function. 


appears. The QSH phase is obtained for all studied films, 
both with odd and also even number of monolayers. How¬ 
ever, for all but the (lll)-oriented SnTe films with an 
even number of monolayers an overlapping of bands in r 
and M diminishes the final band gap. Hence, the edge 
states appear either against the background of the bands 
or within a very small energy gap (see Supplementary 
Section II). The (lll)-oriented SnTe films with an even, 
close to 20, number of monolayers is thus the best candi¬ 
date for observing the QSH effect in the IV-VI TCI. The 
QSH phase produced in this material system is highly 
robust, being even oblivious to the lack of inversion sym¬ 
metry. The latter indicates that a suitable substrate or 
an overlayer will not destroy this 2D TI phase. Quite the 
contrary, using carefully selected substrate or overlay ma¬ 
terial (e.g., Pbi-^Eu^Sei-^Tey quaternary alloy) would 
allow to tune the height of the energy barriers and biax¬ 
ial strain in the structure. 28 As shown already in Ref. 29 
for SnTe (001) surface and for SnTe nanomembrane in 
Ref. 30, the uniaxial or biaxial strain can be used for 
manipulating with the energy gaps in the surface (or in¬ 
terface) states. The strain can be also used for solving 
the mentioned above problem of band overlapping, which 
diminishes the resulting energy gap (as it was shown in 
Ref. 31 for solving a similar obstacle in PbTe/(Pb,Sn)Te 
heterostructures). Thus, these results may pave the road 
for experimental realization of QSH in SnTe and SnSe 
layers and their quantum well structures. 
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